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ABSTRACT 

Long-term X-ray modulations on timescales from tens to hundreds of days have been 
widely studied for X-ray binaries located in the Milky Way and the Magellanic Clouds. 
For other nearby galaxies, only the most luminous X-ray sources can be monitored 
with dedicated observations. We here present the first systematic study of long-term 
X-ray variability of four ultraluminous X-ray sources (ESO 243-49 HLX-1, Holmberg 
IX X-1, M81 X-6, and NGC5408 X-1) monitored with Swift. By using various dy¬ 
namic techniques to analyse their light curves, we find several interesting low-frequency 
quasi-periodicities. Although the periodic signals may not represent any stable orbital 
modulations, these detections reveal that such long-term regular patterns may be re¬ 
lated to superorbital periods and structure of the accretion discs. In particular, we 
show that the outburst recurrence time of ESO 243-49 HLX-1 varies over time and 
suggest that it may not be the orbital period. Instead, it may be due to some kinds 
of precession, and the true binary period is expected to be much shorter. 

Key -words: methods: data analysis - X-rays: binaries - X-rays: individuals 
(ESO 243-49 HLX-1, Holmberg IX X-1, M81 X-6, NGC 5408 X-1) 


1 INTRODUCTION 

Over the last two decades, all-sky X-ray monitoring observa¬ 
tions of X-ray sources in the Milky Way and the Magellanic 
Clouds have revolutionised our understanding of X-ray bi¬ 
naries, leading to many important discoveries. For instance, 
detections of X-ray outbursts of transient sources allow us to 
discover new black hole candidates, and subsequent identi- 
hcation at other wavelengths enable a detailed study of the 
nature of the system. By following the X-ray evolution of 
X-ray transients, one can study different intensity/spectral 
states and investigate the nature of the compact object. In 
addition, some persistent X-ray binaries show X-ray variabil¬ 
ity on timescales from milliseconds to years. While observa¬ 
tions with high timing resolution can probe variability from 
milliseconds to hours, orbital periods and superorbital peri¬ 
ods (periods longer than the orbital periods) on timescales 
of days to years can only be studied with regular monitoring 
observations. 


* E-mail: akong@phys.nthu.edu.tw 


All these observations have been proven to be very suc¬ 
cessful for our understanding of X-ray binaries, but the sam¬ 
ple is mainly limited to the Milky Way which suffers diffi¬ 
culty in determining the distance to the source (hence uncer¬ 
tain luminosit y measurements) and high extinction in some 
regio ns (e.g., iNavlor et al.l Il99ll : [Thompson fc Rothschildl 
I 2 OO 9 II . Apart from a few cases in the Magellanic Clouds, 
monitoring observations of X-ray binaries are very rare in 
other nearby galaxies. With the improvement of the spa¬ 
tial resolution and sensitivity of X-ray telescopes, such ob¬ 
servations have become feasible now. One intriguing dis¬ 
covery from X-ray observations of nearby galaxies is ul¬ 
traluminous X-ray sources (ULXs), which are luminous 
(Lx(o.3-io kev ) > 10^® ergs s“^), non-nuclear X-ray point¬ 
like sources with apparent X-ray luminosities above the Ed¬ 
dington limit for a typical stellar-mass (~ IOMq) black hole. 
The majority of ULXs are believed to be accreting objects 
in binary systems with X-ray flux vari ability on timescales 
of hours to years (e.g, iHeil et aP I 2 OO 9 I for short-term vari¬ 
ability). Assuming an isotropic X-ray emission and that the 
source does not exceed the Eddington limit, a ULX is the 
best candidate for an intermediate-mass black hole (IMBH) 
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with a mass of ~ 10^ — 10 ^Mq (e.g.. iMakishima et alJli^OOd : 
iMiller fc Colber^l2004h . While ULXs may represent a miss¬ 
ing link between stellar-mass black holes and supermassive 
black holes in galactic centers, their formation and evolu¬ 
tion are still not well understood. Unlike Galactic X-ray bi¬ 
naries, ULXs in nearby galaxies cannot be detected with 
typical all-sky X-ray monitor instruments because of their 
X-ray faintness. However, one can observe ULXs with large 
instruments regularly. In the past 10 years, RXTE, Chan¬ 
dra, XMM-Newton and Swift have been used for this pur¬ 
pose. For some well-known ULXs, long-term time series data 
have been building up allowing us to study their orbital and 
superorbital periods. Recently, there are some reports on the 
possible orbital periods of a few ULXs based o n Swift mon¬ 
itorin g observations (e.g., IStrohmaveil l2009l : iLasota et al.l 
l201lh . Due to the limited time span, the stability of the 
period which is one of the important indicators for an or¬ 
bital period requires verification from more observations. 
Moreover, due to a non-regular sampling and instrumental 
artefacts, the results have to be cross-checked with different 
timing analysis techniques. 

In this paper, we report an investigation of long-term 
X-ray variability of four ULXs (ESO 243-49 HLX-1, Holm- 
berg IX X-1, M81 X- 6 , and NGC5408 X-1) that have been 
observed regularly with Swift. Although the long-term light 
curves for some of the selected sources in this paper were 
studied before, we present a systematic timing analysis us¬ 
ing both static and dynamic Fourier power spectra on a 
much longer timeline of data. We also considered the contri¬ 
butions from white noise and red noise to derive the signifi¬ 
cance of a signal. Finally, we applied a wavelet analysis and 
Hilbert-Huang transform to compare with results obtained 
by Fourier analysis. For the present study we concentrate 
on the usage of different timing analysis techniques for the 
study of long-term X-ray variability of ULXs. 


2 OUR SAMPLE 

In this study, we select ULXs {Lx > 10®® erg s“^) that have 
been continuously monitored by Swift for more than 4 years. 
Because Swift is not an all-sky monitoring instrument, the 
sample is limited by the observations in the archive that 
were originally proposed for different purposes. In order to 
perform meaningful timing analyses, we require that each 
target must have more than 200 data points. In total, we 
have four ULXs (ESO 243-49 HLX-1, Holmberg IX X-1, 
M81 X -6 and NGC 5408 X-1) in our sample that satisfy our 
criteria (see Table [T] for an observing log). In the following 
subsections, we briefly introduce our targets. 


Table 1. Summary of observations used in analysis 


Source 

Date of observations 

Data points 

ESO 243-49 HLX-1 

2009 August - 2013 December 

251 

Holmberg IX X-1 

2008 December - 2013 October 

483 

M81 X-6 

2009 April- 2013 October 

396 

NGC 5408 X-1 

2008 April - 2012 August 

371 


(|Lasota et al.l[201ll~) until recent reports for a delay of the 
latest outburst ( Godet et al.ll20i^ : iKong et al.ll2ni5h . 


2.2 Holmberg IX X—1 

Holmberg IX X-1 is a well-known ULX with an X-ray lumi¬ 
nosity of ~ 10"^® ergs s“® in 0.5-10 keV b and. It was firs t 
discovered by the Einstein Observatory llFabbiand Il988l l , 
and has been observed by all major X-ray observa tories 
throug hout the las t 20 y ears IlI Parola et al.l 1200il l. Re¬ 
cently, iKong et al.l (I 2 OI 0 II also used Swift to monitor the 
X-ray evolution of Holmberg IX X-1 and utilised the co¬ 
added spectra taken at different luminosity states to study 
the spectral behaviour of the source. The best spectral fits 
are provided by a dual thermal model with a cool blackbody 
and a warm disc blackbody. This suggests that Holmberg 
IX X-1 may either be a IOMq black hole accreting at seven 
times above the Eddington limit or a IOOMq maximally ro¬ 
tating black hole accreting at the Eddington limit, and we 
are observing both the inner regions of the accretion disc 
and outflows from the compact object. 


2.3 M81 X-6 

M81 X- 6 , which is also known as NGC 3031 X-11, is 
the brightest non-nuclear X-ray source in M81. The aver¬ 
age X-ray luminosit y of the source in 0.3-10 keV band is 
~ 2 X 10®® ergss~^ ilSwartz et al.ll2003ll . It was suggested 
based on optical colours from Hubble Space Telescope obser- 
vations that the system contains a 23Mq 08 V companion 
llLiu et al.ll2002h . Using X-ray s pectroscopy, the bl ack hole 
mass is estimated to be I 8 M 0 (ISwartz et al.l 120031 1. Based 
on the mass of the secondary star, and by assuming that 
it fills the R oche lobe, an or bital period of about 1.8 days 
is estimated llLiu et al.ll2002f l. Although the cadence of the 
current Swift monitoring programme for ULXs is usually 
from a few days to a week and such a short orbital period 
is difficult to be detected, we can test if a longer periodicity 
exists. 


2.1 ESO 243-49 HLX-1 

The ULX ESO 243-49 HLX-1 is the most luminous ULX 
reaching a peak X-ray luminosities of 1.1 x 10“*® ergs s“® 
in 0.2 - 10 keV band dFarrell et al.l [2009h . Assuming an 
isotropic emission in a super-Eddington state, a conserva¬ 
tive minimum black hole mass is estimated to be ~ 5 OOM 0 , 
making it a very promising IMBH candidate. During the 
past few years, the source exhibits several X-ray outbursts. 
The regular separation between each outburst of ESO 243- 
49 HLX-1 is assumed to be caused by the orbital motion 


2.4 NGC 5408 X-1 

NGC 5408 X-1 is one of the most luminous ULXs with an X- 
ray luminosity of ~ 2 x 10“^® ergs s“® in 0.3-10 keV band. It is 
also one of the hrst ULXs observed regularly with Swift. Us¬ 
ing t he Swift/X-'Ray Te lescope (XRT) data spanning ~ 500 
days. I^rohmaverl (l2009ll discovered a 115.5±4-day periodic¬ 
ity and suggested that it is due to orbital modulation. If this 
interpretation is correct, it would support that NGC 5408 X- 
1 contains a ~ IOOOM 0 black hole wi th a 3 — 5 M 0 compan- 
ion based on a theoretical simulation (IPorteeies Zwart et al.l 
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Figure 1. Swift XRT 0.3-10 keV light curves of ESO 243-49 HLX-1, Holmberg IX X-1, M81 X-6, and NGC 5408 X—1 used in our 
analysis. Data points before the determined starting epoch were eliminated to avoid any fake signals arising from huge data gaps. The 
red dashed and blue dotted lines are used to label the interval of those intriguing detections summarized in Table [2] 


I 2 OO 4 I I . iHan et 51 ll2012ll also used the weighted wavelet z- 
transform to confirm the periodicity of 115.5 days, but they 
only cons idered the li g ht cu rve within the first 500 days. 
However, [Foster et al.l (120101 1 argued t hat such a modula - 
tion is due to a superorbital period Isee lCharles et ahlbOOSl l 
by comparing its physical properties to the micro-quasar, 
SS 433. They suggested that the 115.5-day modulation is 
originated from a precessing jet, and this scenario implies 
that the system has a much shorter orbital period and 
the accreting object is a stellar- mass black hole. Using 4 
years of observations from Swift dCrise et al ] |2013h . an or- 
bital origin of the ^ 115-day periodicity can be rejected. 
IPasham fc Strohmaveil (l2013l l revisited their interpretation 
as a quasi-sinusoidal X-ray modulation, and re-determined 
an orbital period of 230 days for NGC 5408 X-1 based on 
the average time scale for the recurrence of sharp dips in the 
X-ray intensity. 


3 DATA REDUCTION AND ANALYSIS 
3.1 Data Reduction 

All the data reported in this paper were taken with the X- 
ray telescope (XRT) onboard Swift. We only used data taken 
in the photon counting mode so that we can identify the pre¬ 
cise locations of our targets without any contamination. We 
extracted XRT light curves by usi ng the XRT products gen¬ 
erator Q (lEvans et al.ll20()7l . l2009ll . It first creates an image 
from the event list and identifies our target for each obser¬ 
vation. Only events with energy in the range 0.3-10 keV 
with grades 0-12 are included. A circular source extraction 
region is chosen to match the point-spread-function. For the 
background, an annulus centered on the source is used. Then 
for each observation, source and background counts are ex¬ 
tracted. Source count rates are corrected for the good time 
interval, and losses due to bad pixels and bad columns. 

Scatter data points can be present at the beginning of 
each Swift monitoring. These data are not taken into consid¬ 
eration in our analysis, in order to avoid spurious effects due 


^ http://www.swift.ac.uk/user_objects 
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to large observational gaps. We also applied a criterion to 
reject data with nncertainties of the background larger than 
3 (t. In addition, we removed those data points with anoma¬ 
lously high or low count rates to avoid contamination by 
the high background or spurious flaring events. Some data 
points with count rate errors larger than their count rate due 
to very short exposures 10 — 30 s) were not included in 
the temporal analysis as well. Fig. [T]shows the original light 
curves of all four ULXs in our sample. Since the window 
function may generate artificial signals, we also re-binned 
the data samples every 5 days in our analysis to compare 
with those detections obtained from the original light curve. 
Because the mean cadence of observations is about 2~10 days 
(excluding obvious data gaps), we only considered periodic¬ 
ities longer than 10 days (i.e., < 0.1 1/day in frequency) in 
our timing analyses. 


3.2 Methods 

In order to look for any modulation, we first performed 
a timi ng analysis by using the Lomb-Scargle periodogram 
fLSP: II^mblll97^ : IScargIg 1982l ~). a modification of the dis¬ 
crete Fourier transform which is generalized to the case 
of unevenly spacing. Our static periodogram/power spec¬ 
tra were generated through the REDFIflo software and cor¬ 
rected for the bias arising from unevenly spaced data in the 
time domain that overestimates the high-frequency power 
llSchulz &L Statteggeilll997l : ISchulz k. Mudelseell2002 i. Once 
period candidates are found, we examine their long-term sta¬ 
bility by using a dynamic power spectrum (DPS). Finally we 
cross check the results with a weighted wavelet ^-transform 
(WWZ) that has been known as an excellent indicator of 
signals at the expense of being less sensitive to their powers. 

In the following, we will explain the detailed setting 
of our methods. To derive the LSP using REDFIT, the 
oversampling (OFAC) and high-frequency (HIFAC) factors, 
which set the period range and resolution ( Press fc Rvbickil 
Il989h . were individually set as 64 and 2 for each process to 
cover the frequency range in our analysis. In order to com¬ 
pare with the power determined by a typical LSP, we trans¬ 
lated the obtained amplitude of the spectrum into Lomb- 
Scargle power dividing by a scaling factor related to the 
average sampling rate and the variance of count rate in 
our data sets. The false alarm level of the white noise was 
set equal to 1 — (1 — where P is the highest peak 

in a periodogram and Ni is the number of independent 
trial frequencies. In order to simplify our computation, we 
adopted the empirical fu nction obtained from simulations 
dHorne fc Baliunaa Il986| l to describe th e number of inde - 
pendent frequencies (Horne number; e.g.. lKong et al.lll998ll . 
We also fit a first-order autoregressive (ARl) process to the 
time se ries to estimate the red n oise spectrum (solid line in 
Fig. (2] ISchulz fc Mudelseel I 2 OO jl . The false alarm level of 
the red noise was assessed through scaling the theoretical 
red-noise spectrum by an appropriate percentile of the x^~ 
probability distribution. The analysis with LSP has been 
widely applied for studying long-term X-ray variability in 
X-ray binaries fe.g.. lSmale fc Lochnerlfl992l : IWiinands et al.l 


^ http://www.ncdc.noaa.gov/paleo/softlib/redfit/redfit.html 


1 19961 : iKong et al.lll998l : iKotze fc Charles! l2012fl . and RED- 
FIT was also successfully applied to the power spectrum to 
investigate the red- noise level of long-term X-ray variability 
dFarrell et al.|[2009h . 

For the DPS , we em ployed a similar method outlined 
in IClarkson et ahl (l2003al lal. The obtained archival data of 
our targets were analysed with a sliding window to generate 
the power density, and this window was moved to cover the 
entire time series. This method accounts for variations in 
the number of data points per interval, while instability of 
quasi-periodicities can be resolved with a DPS. The size of 
the window depends on the length of the periodicity to be 
examined in the DPS. The minimum size of the window was 
set such that it covers at least one complete cycle. Since 
suspected periodicities of our detections range from ~ 30 to 
300 days, the window length was set to 500 days. We set the 
step size of the sliding window in the DPS at ~ 10 days to 
avoid artificial effects. 

The code implementation of WWZ that we used is 
based on a theoretical study bv lFoste'3 (ll996cl L For an ob¬ 
served time series with N data sets (e.g., x{ta)) taken at 
N discrete times {ta : a = 1,2,..., N}, we can define N- 
dimensional contravariant vector in sampling space with a 
representation of the canonical basis. Any function of time 
f{t) can also be defined as a contravariant vector in sam¬ 
pling space, and therefore the inner product of two functions 
f{t) and g{t) is defined as 


<f\g>= 


Wc/(tc)g(tc) 

E N ’ 


( 1 ) 


where Wa is the statistical weight assigned to uneven data 
points. For a set of trial functions < : a = 1, 2,..., r >, 

the S'-matrix for the sampling space is given by 


Sab =< > 


( 2 ) 


and the coefficients of equally weighted data can be de¬ 
scribed as 

= (3) 

b 

for the model function. 

The wavelet transform (WT) onto the complex trial 
function at the frequency u: and the time r is 

/(i) = (4) 


where c is a constant, and can be regarded as a projection 
onto the complex trial function of with real sta¬ 
tistical weights chosen as Wa = Ua-T) WWZ 

consists of the effective number and the weighted variations 
of uneven dat a and model functions as below. According to 
iFosteil lll996al|H) , the effective number is 








( 5 ) 


The weighted variation of the uneven data x is given as 


K = 


Ea '^olX (ta) 

T.X'^C. 


Ec 'U)aX{tcf 


Ea^c, 


=< x\x > — < l\x > 


and the weighted variation of the model function y is 


Ea^^^ 


Ec«'^“g(tc«) 


Ea^e, 


=< y\y> - < i|g > 


( 6 ) 

( 7 ) 
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Then the WWZ is read as 


Z = Z{LO,r) 


{N,ff-3)Vy 
2(V - Vy) 


(8) 


Notably, the coefficients Pa are directly computed based on 
the uneven-smapling data. In practice, the choice of 
could be 


$ 1 ( 1 ) = 1, ^ 2 (t) = cos(oj(t — t)), 'l? 3 (t) = sin(a;(t — r)) (9) 

Once a period candidate is determined from the WWZ, 
the amplitu de of the sign al can be determined as WWA = 
'/^a =2 Va itFosteririOOBcJ l. All the obtained signals in this 
paper are subtracted by their mean, and the significance 
can also be examined with WWA or the corresponding ef¬ 
fective number (eq.[5]). Here we concentrated on the signals 
resolved from the evenly sampled data by the WWZ in a 
comparison to those obtained by the DPS. Comparing with 
other time-frequency analysis methods, the WWZ is opti¬ 
mised for non-stationary unevenly time-series analysis, so 
this method is quite adequate to verify any signal that ex¬ 
ists in the obtained data sets. 


4 RESULTS 

4.1 ESO 243-49 HLX-1 


The long-term light curve of the source is shown in Fig. [T] 
It is clear that the source exhibits substantial variability. In 
particular, there are nearly hve complete outbursts with a 
fast-rise exponential-decay profile similar to typical Galac¬ 
tic X-ray transients. The Brst and second luminosity peaks 
for ESO 243-49 HLX-1 were in A ugust o f 2009 and late- 
August of 2010 llGodet et al.ll20l3 : lKonsll2010l '). With the 
regular monitoring by Swift, the next two cycles of outbursts 
recurred in 2011 and 2012 (iKondlioill : iGodet et af]|201ll : 
iKondlioi^ l. The separation between outburst peaks changes 
from ~ 350 days for the first two outbursts to ~ 370 — 380 
days for the following ones. The most recent outburst was 
detected in early-October of 2013 wit h a delay by almos t a 
month compared to the previous one llGodet et al.ll2oi^ . 

By performing a LSP analysis, the strongest peak higher 
than both 99% red and white noise significance levels for 
binned data is 378 ± 6 days (Fig. [2]), where the error was 
estimated by propagating two uncertainties. The first un¬ 
certainty in frequency of the trial periodicity was estimated 
using 


Sf 


3 1 

8tVp’ 


( 10 ) 


where T is the tota l time interval and P is the peak power 
dLevine et al. I [m3). The second uncertainty accounts for 
measurement errors effects, and was determined by carry¬ 
ing out Monte Carlo simulations. We simulated 10^ light 
curves by randomizing the count rate according to a Gaus¬ 
sian distribution with standard deviation equal to the count 
rate uncertainties. Then, we repeated the LSP analysis and 
derived a standard deviation for the peak frequency of the 
periodicity. The major signal also has a power higher than 
95% confidence level determined from the red noise model 
for the power spectrum generated from the original data, in¬ 
dicating that this period is statistically signihcant. Similarly, 
the DPS in Fig. [3] obtained from the re-binned light curve 


shows a significant periodicity at about 380 days. There are 
two other peaks (190±3 days and 126.2±1.1 days) in Fig. [2] 
with powers more significant than the 99% white noise con¬ 
fidence level, and they might correspond to different har¬ 
monics of the strongest signal as we will descr ibe in detai l 
in section §5.1. According to previous reports (iKon a [mil; 
iGodet et al.ll201ll : lKong|2012l L the separation between the 
peak luminosity of the first two outbursts and the following 
outbursts has increased. Except for an obvious delay of the 
outburst occurred in late 2013, because the poor resolution 
provided by the DPS cannot discriminate such a change un¬ 
ambiguously, we considered the LSP with the sub-divided 
data sets of day < 865 (the black LSP) and day 365-1210 
(the red LSP) such that both include three outbursts (see 
Fig. 2]). The strongest signals shown in Fig. 0] are 341 ± 10 
days and 378±13 days. The difference can be associated with 
the change in the separation between the two outbursts. In 
addition, we found that the power of the signal at ~ 190 
days significantly increases in recent observations. This re¬ 
sult can directly be linked to the second harmonic/first over¬ 
tone with the separation between recent outbursts (~ 380 
days) and provides evidence for variation of the duty cycle 
for each outburst. This can clearly be seen in Fig.[3]as well, 
and a bifurcation of the main periodicity ~ 380 days can be 
found after day ~ 815. This might also cause the delay on 
the most recent outburst in 2013 early-October. Because of 
the lack of observations between day 1210-1435, our results 
were seriously smeared in this time interval. 

The aforementioned phenomena discovered in the DPS 
can also be seen from the WWZ periodogram (see Fig. 0- 
There are three strong signals, which represent periodicities 
of ~ 380 days, ~ 190 days and ~ 126 days. The major 
signal (~ 380 days) across the whole data is very strong, 
but the resolution of the WWZ periodogram is not good 
enough to resolve any variations. However, we can obviously 
find that two sub-signals generated from the overtones only 
get significant after day 720. This result totally agrees with 
those from the DPS, indicating a variation of the duty cycle 
of each outburst. 

Since ESO 243-49 HLX-1 has a regular monitoring be¬ 
fore day 1210, it is an ideal object to be further studied 
with the Hilbert-Huang transform fHHT: [Huang et al.lll998h 
since HHT analysis requires an evenly sampled data set. 
Though the data were binned with 5 days per bin, there 
is also an obvious gap after the 4th outburst between day 
1210-1435 (only two data points were recorded on days 
~ 1353 and 1409). Because this gap is in the quiescent state 
and the count rate is expected to be less than 0.003 cts/s 
according to previous observations, we re-interpolated the 
light curve usin g a piecewise cubic Hermite spline function 
(lKrevszisll2005l ) such that the data are evenly sampled with¬ 
out any significant artificial effect. Following the similar pro- 
cess to study th e superorbital period of SMC X-1 with HHT 
(IHu et al.l 1 20 ni l. we ap plied an ensemble empirical mode 
decomposition (EEMD; IWu fc Huanel[2009l l method to de¬ 
compose the original light curve into a number of intrinsic 
mode functions (IMFs) with the medium noise level at first. 
The decomposed component of the largest energy/power 
contains four cycles, which correspond to the main modula¬ 
tion obtained from the LSP and DPS. Here, w e employed a 
normalised HHT fFig.l^ [Huang fc Londl^OOSlL which is sen¬ 
sitive to both inter-wave modulation (i.e., the variation of 
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ESO 243-49 Holmberg IX X-1 



Frequency (1/day) Frequency (1/day) 



Frequency (1/day) 


Frequency (1/day) 


Figure 2. Lomb-Scargle periodograms of the four ULXs. The black and the red curves represent the bias-corrected LSPs of the unbinned 
data and the 5-day average data, respectively. The red noise models for each power spectrum obtained from the unbinned and re-binned 
light curves are individually plotted as black and red solid lines for comparison. The dashed horizontal and dashed-dotted lines in each 
panel indicate the 99% white and red noise significance levels for the unbinned light curves (in black) and the re-binned light curves (in 
red). 


a cycle length) and intra-wave modulation (i.e., the instan¬ 
taneous frequency changes within one oscillation cycle), to 
represent the instantaneous frequency on the Hilbert spec¬ 
trum and to compare with results obtained from the DPS 
and WWZ. A strong variable signal between ~ 250 and 500 
days (frequency ~ 0.002 — 0.004 1/day) is seen in the Hilbert 
spectrum, and it fits well to the major signal yielded from 
the DPS. The variation shows roughly five cycles, which is 
consistent with the quasi-periodicity detected from the light 
curve. Since the ~ 380-day modulation is highly non-linear 
(non-sinusoidal), the intra-wave modulation would dominate 
the variation of instantaneous frequency. 

4.2 Holmberg IX X—1 

Fig. [Hand Fig.[2]show the long-term light curve and the LSP 
of Hol mberg IX X-1, respectively. Comparing to lKong et al.l 
(l2010l) . Holmberg IX X-1 has gone through at least three 
other high/low intensity tra nsitions. The two peri od candi¬ 
dates (~ 28.8 and 64.5 days; iKaaret fc Fenal[2009l ) found in 
previous analyses below the 85% white/red noise confidence 


level in the LSP. Instead, we obtained a significant peri¬ 
odic signal from both unbinned and re-binned data sets at 
625 ± 20 days, which reaches the 99% red noise significance 
and exceeds over the 99% white noise significance levels in 
Fig. [2] According to the light curve, the peak (~ 0.4 cts/s) 
recurs at days 440-455, 1120-1130 and 1730-1740, and each 
separation is about 625 days. It can also be seen at days 
~ 645 and 1270-1275 with a separation of ~ 625 days. If we 
consider the dip features determined from the light curve, 
we can also find them with a separation of about 600 days at 
days 135-140, 750-760, and 1325-1330. However, the length 
of the current light curve used in our analysis is ~ 1770 
days, which only includes two complete cycles for a period¬ 
icity of ~ 625 days. Given the scattering of the data and 
the total length of the observations, the current data set 
is not suitable for the HHT. Since there is an obvious gap 
without any data points at days 845-1085, we divided the 
data into two groups to inspect whether there is any periodic 
signal in a specific time interval. These results are also pre¬ 
sented in Fig.m Because the WWZ can also be employed to 
unevenly sampling data, we double checked those signals re- 
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Time (Day since 2009 August 5) 



Time (Day since 2008 December 1) 



Time (Day since 2009 April 1) 



Figure 3. Dynamic power spectra of ESO 243-49 HLX-1 (upper left), Holmberg IX X-1 (upper right), M81 X-6 (lower left), and 
NGC 5408 X-1 (lower right). The DPSs were obtained with dynamic Lomb-Scargle method with a moving window of 500 days using the 
re-binned data. The colour of the contour represents the strength of the LSP power. 


solved from the DPS and from the LSP with the sub-divided 
data. In the WWZ periodogram of Fig. [6] there is a strong 
signal at low frequencies (~ 625 days) persisting across the 
whole data set. The significance of this signal slightly de¬ 
creases after day 1050, and a similar phenomenon can also 
be observed in Fig. [3] We cannot confirm such a long-term 
periodicity with the current data with limited cycles since 
the resolution of signals obtained from the WWZ is not good 
enough. Furthermore, no adequate conclusion can be made 
on this feature from DPS because the adopted window time 
(500 days) is shorter than the periodicity we are discussing. 

4.3 M81 X-6 

We show the Swift long-term light curve of M81 X-6 in 
Fig. m The LSPs obtained from both the unbinned data 
and 5-day average data are presented in Fig. [2] The source 
varies by a factor of ~ 4 on a time scale of months. The 
LSP indicates that the strongest candidate periods are at 
about 110 days and 370 days, but they are not statistically 
significant (< 99% red noise confidence level) when we take 
into account the entire data sets. 

Only two signals are more significant than the 99% 
white noise confidence level in the DPS (Fig. [3]), and this 
is consistent with the LSP in Fig. [T] The black curve de¬ 
scribes the LSP obtained from the re-binned data sets of 
day 740 and the red one represents the LSP obtained 


from the re-binned data sets of ~ day 965-1645. The pe¬ 
riod at ~ 370 days is only prominent before day 740 while 
the period at ~ 110 days is only significant after day 965. 
There are two obvious gaps without any data points in the 
light curve after day 965. The length of the two gaps (at 
days 1160-1315 and 1350-1490) is 155 days and 140 days. 
In addition, we also found that the duration with frequent 
cadence at days 965-1160 and 1490-1645 is 195 days and 
155 days, respectively. The duration of the gaps and the 
data with frequent cadence also propagate the fake signals 
in LSP of Fig. (2] When we examined the signals with the 
WWZ shown in Fig. [S] the distinction between the two ma¬ 
jor signals and others are much more evident. These findings 
strongly suggest that those sub-signals found in the LSP of 
Fig. [5] are artificial effects. 

4.4 NGC 5408 X-1 

With Swift/yiRT data covering about 1500 days, we plot the 
long-term light curve in Fig[T] We performed a power spec¬ 
tral analysis with LSP and the significance levels of white 
and red noises in the resulting power spectrum are shown 
in Fig. [2] We do not find any significant periodicity; the 
highest peaks obtained from the original data and the re¬ 
binned data sets (5 days per bin) are at 112.1 ± 1.0 days and 
189 ± 3 days, respectively. These results are consistent with 
the proposed (112.6 ±4)-day quasi-periodic modulation ob- 


© 0000 RAS, MNRAS OOO.ITlITSl 











8 L. C.-C. Lin, C.-P. Hu, A. K. H. Kong, D. 


C. C. Yen, J. Takata and Y. Chou 



Frequency (1/day) 



0 I V I n t V 

0.002 0.005 0.01 0.02 0.05 0.1 


Frequency (1/day) 




Frequency (1/day) 


Figure 4. Lomb-Scargle periodograms of the four ULXs at different epochs (see Section 4.1-4.4 for the definitions of different epochs). 
Re-binned light curves were used to derive the LSPs. The power spectrum of each epoch is marked with different colours, while the 
dashed straight lines and dashed-dotted lines represent the corresponding 99% white and red noise confidence levels of the LSP power 
at each epoch. Details of the adopted time intervals are described in the main text. 


tained bv IPasham fc St rohmay^ (|2013l) and the ~ 187-day 
periodicity detected by Grise et al. (1201^) . Because there is 
a large deviation between the obtained peak at ~ 190 days 
and the proposed orbital period ( 243 ± 23 days) of this sys¬ 
tem JPasham fc Strohmaveiir2013l ). it is difficult to establish 
a direct connection between our results and the suggested 
orbital modulation. 


In the DPS (Fig. [S]), it is clear that the signal near 112 
days is quite strong during the first ~ 500 days of obser¬ 
vations, but it has become weaker after that. The quasi- 
periodic signal near 112 days has significantly decreased af¬ 
ter day 800 and evolved into a quasi-periodicity of ~ 190 
days. On the other hand, the ~ 190-day period has been 
evident after day 800 and was quite weak before day 500. In 
order to clarify this, we used the LSP obtained in different 
time intervals to determine the strength of signals obtained 
in different stages and to avoid any contaminations from ar¬ 
tificial signals arisen from the length of empty observations. 
The black and red curves in Fig. [4] represent the LSP re¬ 
sulted from the re-binned data of ~ day 0-535 and ~ day 
0-1230 since MJD 54565.8, respectively. An exchange of the 
significance between two major signals (~ 112 days and 190 
days) is clearly shown although they are not statistically sig¬ 


nificant, and the similar phenomenon can also be verified by 
the WWZ shown in Fig. [6] 


5 DISCUSSION 

We generated LSP and DPS for four ULXs from the long¬ 
term monitoring Swift XRT data, with the aim of revealing 
periodic signals embedded in the light curves. The results 
were also examined with the WWZ. In the case of ESO 243- 
49 HLX-1, we also used HHT to further study its periodicity. 
Here, we discuss the results for each object. 


5.1 ESO 243-49 HLX-1 

Both the LSP and DPS shown in Fig. [2] and Fig. [3] present 
a strong persistent periodic signal at ~ 380 days. This de¬ 
tection is consistent with the mean separation between the 
start of different outbursts or between their peak luminosi¬ 
ties (~ 350, 373, 375 and 407 days). The gradual delay of 
each outburst also reflects on the result of DPS, and a clear 
bifurcation from the main signal can also be seen. Because 
the biforked feature is only significant in the data > day 
720, we divided the data sets observed with similar cadence 
(<~ 1210 days) in two segments that both include the entire 
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Figure 5. Comparison between DPS and HHT for ES0 243—49 HLX-1. The DPS obtained with a moving window of 500 days from the 
re-binned light curve (contours) is overlapped on the Hilbert spectrum. It is clear that they are consistent with each other. 


duty cycle of three outbursts to trace the variation of the 
periodic signals. The LSP in Fig. 3] yielded from the data 
which include outbursts from the second to the fourth re¬ 
veals some other quasi-periodic signals that correspond to 
the harmonics of the main signal at ~ 380 days. If we re¬ 
duce the length of the light curve to include the third and 
the fourth outbursts only, we find that the significance of 
these signals is enhanced. 

The Hilbert spectrum and the WWZ provide us an op¬ 
portunity to verify the results yielded from LSP and DPS. 
We compared the DPS and the Hilbert spectrum in Fig. 
The contour represents the power of DPS generated with 
a moving window of 500 days through the re-binned light 
curve. The colour map represents the Hilbert energy de¬ 
fined as the square of the amplitude, and it was smoothed 
using a Gaussian filter. It is clear that the Hilbert spectrum 
yielded from the normalised HHT passes through the peak of 
DPS, indicating that the result of the Hilbert spectrum and 
the DPS are consistent. Furthermore, the Hilbert spectrum 
provides a much higher/better resolution to investigate the 
structures of the variability/instability caused by the intra¬ 
wave modulation, i.e. the frequency modulation within one 
cycle. It is caused by the non-sinusoidal nature of the mod¬ 
ulation because the modulation shape contains a sharp rise 
and an exponential decay. The maximum frequency of the 
the main signal due to the intra-wave modulation appears 
at ~ 390, 745 and 1125 days (Fig. [5]), which are similar to 
the epochs of the peak luminosity of each outburst. 

Comparing to the DPS in Fig. [3] we do not find a bi- 
forked feature on the WWZ periodogram. The reason is that 


the Lomb-Scargle method tends to describe the data with 
a smooth sinusoidal wave characterized by a longer period 
when the moving window approaches the flat low count rate 
state. Nevertheless, the gradual delay of the major signal can 
still be seen in Fig. [S] if we zoom in the frequency domain. 
We can also find that a new signal begins to appear after 
day 510, and its power is enhanced after day 760. Because 
day 510 is at the quiescent state after the second outburst 
and day 760 is close to the peak of the third outburst, we 
speculate that the emergence and the enhancement of this 
new signal might be due to a change of the shape of the 
light curve in outbursts. The best fit to the cycle length of 
the first four complete outbursts can be estimated as ~ 370 
days from the strongest signal in the LSP. We then folded 
the light curve with this periodicity to examine the differ¬ 
ence of the outburst light curve profiles in each cycle. As 
shown in the top panel of Fig. [7] we use different colours 
to mark each cycle. The black curve characterizes the first 
cycle, while the blue one represents the second cycle. The 
third and the fourth outbursts are depicted with green and 
red colours, repectively. The most recent outburst denoted 
by the yellow curve recurred with a relatively long separa¬ 
tion 407 days) and clearly shows a delay during its rising 
phase. The first four outbursts start with a sharp rise near 
phase 0, but the profiles of their decays are different. The 
duty cycle of the outburst decreases from ~ 50% in first 
cycle to ~ 30% in the last cycle. The third and the fourth 
outbursts have similar phase for each peak, and a good fit to 
their duty cycle can be modelled by a consecutive sinusoidal 
signal with a periodicity of 190 days as shown in the bottom 
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Figure 6. WWZ periodograms of ESO 243-49 HLX-1 (upper left), Holmberg IX X-1 (upper right), M81 X-6 (lower left), and NGC 5408 
X-1 (lower right). All the periodograms were obtained with the WWZ stated in §3.2 using the re-binned data. The colour of the contour 
represents the strength of the significance defined through eq. |8] 


panel of Fig. [7] This periodicity is consistent with the signal 
(190 ± 3 days) that we obtained in Fig. (2] 

A third signal of ~ 125 days in Fig.[4]is roughly consis¬ 
tent with the length of the duty cycle of recent outbursts. 
These results give an interpretation to the bifurcation of 
the main signal shown in the DPS and the emergence of all 
the sub-signals shown in the DPS and the WWZ, and these 
quasi-periodicities are indeed originated from the harmonics 
of the main signal. Because the distribution of the data is 
highly non-linear and the duty cycle of the third and fourth 
outburst is much shorter, the EEMD tends to decompose the 
peak/amplitude into another IMF and this feature can also 
be seen in the Hilbert spectrum. Therefore the power/energy 
of the main signal during the Erst two outbursts shown in 
the DPS is weaker than that during the next two outbursts 
(Fig. [31) while the Hilbert spectrum presents an opposite re¬ 
sult. Moreover, we see an obvious signal varying between 
125 and 250 days at day 510-1310 in the Hilbert spectrum 
as well. 

Since the main periodicity embedded in the light curve 
of ESO 243-49 HLX-1 arises from the separation of each out¬ 
burst, this signal is not stable because of the delay between 


outbursts. We can clearly see this feature in Fig. |31 and the 
main period increases to ~ 380 days if we take into account 
the latest outburst. To explain the ~ 380-day modulation, it 
has been suggested that the companion s tar is in a highly ec¬ 
centric orb it with the central black hole (iLasota et al.|[2()lll : 
ISoriallioi^ ). Our analyses, however ^ indicate that the mod- 
ulation is not stable. More recently. lOodet et al.l (120141 ) pro¬ 
posed that the orbital period could change between difier- 
ent orbits as a consequence of stochastic fluctuations. As 
an alternative explanation we speculate that this modula¬ 
tion represents a superorbital periodicity, while the orbital 
perio d of the system i s likely < 125 days (iKotze fc Charle^ 
I 2 OI 2 I) . [King fc Lasotal (|2014l ) suggest that the system may 
be similar to the Galactic microquasar SS433 for which the 
X-ray modulation is due to precession of an X-ray beam. 

5.2 Holmberg IX X—1 

iKaaret fc Fend ll2009h detected two strong possible signals 
near 28.8 and 64.5 days with the Swift data taken between 
2006 July and 2009 June. However, both signals are not sta¬ 
tistically significant when the 11 observations before MJD 
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Figure 7. Folded light curve and sinusoidal fit to the sub-signal obtained from ESO 243—49 HLX-1. Top: The unbinned light curve was 
folded with 370 days to demonstrate the recurrence of each outburst. The phase zero was set according to the first data point used in 
our timing analysis. This light curve includes 5 outbursts recorded by Swift. The recent outburst has a significant delay, which can be 
clearly seen from the rising phase of main peak. Bottom: The light curve taken from the third to the last outbursts overlapped with a 
sub-signal of 190 days. This signal provides a good fit to the duty cycle of the third and the fourth outbursts as shown in Fig. [2] and 

Fig.H 


54200 are removed and only 177 observations after MJD 
54801 (day 878 in Fig. [T]) are taken into the consideration 
dKong et al.ll2O10h . Here, we have more observations cover¬ 
ing ~ 1770 days to investigate the long-term periodic signal 
of Holmberg IX X-1. By taking red noise into account, the 
two periods are not significant. 

We also find a prominent signal corresponding to a 
quasi-periodicity of ~ 625 days (~ 0.016 1/day in Fig. [2]). 
We find that this signal roughly corresponds to the separa¬ 
tion between each peak or dip event although some of them 


seem to appear randomly. According to Fig.|3]and Fig. [4] the 
signal of ~ 625 days is much more prominent in the begin¬ 
ning (before day ~ 845) of the re-binned light curve, but the 
length of this segment is shorter than two complete cycles. 
Although this long-term signal has a stronger significance in 
the WWZ (see Fig. [6ll, the length of the entire light curve 
we used in the analysis is too short (only about 2.8 cycles) 
to confirm this periodicity. In addition, we must emphasize 
that the window size adopted in the current DPS is only 500 
days and it is also the reason why we only show the signals 
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Phase 



Superorbital Phase 

Figure 8. Folded light curve obtained from the Swift data of M81 X-6 between day ~ 965 — 1645. Top: The unbinned light curve was 
folded with 112.4 days. The phase zero (at ~ day 966.5) was set at the epoch of the first observation in the second segment of the light 
curve. Different colours are used to show the data of each cycle. We ignored the 3rd cycle without any observation and the final cycle 
with only 2 data points. Bottom: Binned light curve of the last segment folded with the 112.4-day quasi-periodic signal. 


of interest with periodicity < 500 days in most LSPs. All the 
above reasons make the confirmation of such a long period¬ 
icity difficult. Because the light curve of Holmberg IX X-1 
is quite spiky, it is not appropriate to interpolate the data 
sets to generate an evenly sampled light curve for HHT. In 
order to confirm the long-term periodicity of Holmberg IX 
X-1, more observations with similar cadence are required in 
the future. 


5.3 M81 X-6 

For M81 X-6, except for a major gap with only two observa¬ 
tions at day ~ 740 — 965 and two other gaps with no obser¬ 
vations at ~ 1160 — 1315 and 1350-1490, the mean cadence 
of the data set is 2.8 days. Therefore, the Swift light curve 
is not ideal to e xamine the proposed 1.8-day orbital period 
llLiu et al.ll2002l) . but it provides us an opportunity to study 
the long-term variability of this high-mass X-ray binary sys¬ 
tem. Although there is no significant period in the LSP 
(Fig. [2]), there are two significant periods when we consider 
the light curves in different epochs (Fig. |4]). From the DPS 
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(Fig-O and the LSP obtained from different stages shown in 
Fig. [4] we find that the ~ 370-day period is only prominent 
in the first segment (<~ day 740) of the light curve while 
the ~ 110-day one is only prominent in the last segment (~ 
day 965-1645) of the light curve. The WWZ periodogram 
also clearly demonstrates two strong quasi-periodicities ex¬ 
isted in a different time interval. Comparing with the DPS, 
WWZ is unaffected by the noise generated from the window 
effect and all spurious sub-signals are less significant in the 
WWZ map. 

In the first segment of the light curve, the dip events 
(< 0.015 cts/s) take place at days ~ 170 and ~ 555 and the 
peak events (> 0.078 cts/s) appear at days ~ 25 and ~ 380. 
The separation between the dip or the peak luminosity in 
the light curve provides a possible explanation to the quasi- 
periodic signal at ~ 370 days. Because the length of the 
light curve in the first segment is only about 740 days, which 
only covers about two cycles of the ~ 370-day quasi-periodic 
signal, it is difficult to confirm the stability of this signal with 
this limited data set. 

The ~110-day quasi-periodic signal is not only the most 
significant one we obtained in the second segment of the light 
curve, but also the most significant one in the unbinned 
data. The light curve folded with ~ 110 days is shown in 
Fig. [S] Although there are two obvious gaps in the 3rd cycle 
and between the 4th and the 5th cycles, observations in the 
last segment folded with the ~ 110-day quasi-periodic signal 
show similar structure as demonstrated in the binned folded 
light curve. According to the DPS and the WWZ, this sig¬ 
nal is still significant from more recent observations. Future 
observations with similar cadence are required to infer the 
nature of this quasi-periodicity. 


5.4 NGC 5408 X-1 

With a long-term monitoring of NGC 5408 X-1, we find that 
the most obvious signals from LSP and DPS are ~ 112 days 
and ~ 190 days. However, none of them are statistically 
signifi cant, and this result is consistent with iGrise et al.l 
(l2013tl . In the analysis with DPS, the ~ 112-day period is 
stronger during the first ~ 500 days of observations, which 
is consistent with the (115.5 ± 4)-day period detected by 
IStrohmaver fc Mushotzl^ (120091 ') using the data of the first 
~ 485 days. Another quasi-periodicity of ~ 190-day starts 
to appear after day 400 and becomes prominent after day 
800. 

We do not find any strong signal with a period longer 
than 200 day s to support a proposed or bital period of 
243 ± 23 days (IPasharn fc Strohmaverj[2013f ). Another weak 
signal of ~ 310 days arisen at day 600 shown in Fig. [2] and 
Fig. [^is not related to the proposed orbital period as well, 
even if we take into account its 95% uncertainty. Since the 
recurrence of the X- ray dips related to the claime d period¬ 
icity of ~ 243 days (IPasham fc Strohmaverj|2013f) only oc¬ 
curs before day 1200 and the time duration of each cycle is 
quite unstable (251 days, 211 days, 276 days and 235 days), 
this signal is unlikely real and future monitoring observa - 
tions may see a suppression of the signal. iGrise et al.l (l2013l f 
also indicated that the dip events do not repeat at every 
cycle. Such phenomenon can also be seen in some Galac¬ 
tic tr ansient sources (e.g., GRO J1655-40; iKuulkers et al.l 
I2OO0I I. Because the transient behaviour can cause a change 


on the structure of the accretion disc, the efficiency in X- 
ray irradiation of the disc will change ac cordingly caus¬ 
ing d ifferent behaviours of dipping activities (IKuulkers et al.l 
l2000lf . However, NGG5408 X-1 is a persistent source with¬ 
out any obvious change in X-rays on timescale of years. 
A lack of dipping activity in some epochs has been seen 
in Galactic X-ra y dipping sources (e.g., X1916-053 and 


XB 1254-690, see 


ISmale fc Wachteil 


Chou et ^ I2001I : iHomer et al.l I2001I and 


I999I ) and it may be due to a decrease of 


the vertical structure of the accretion disc. Alternatively, a 
tilted and precessing disc may provide a possible scenario on 
the v ariability of the period icities and the disappearance of 
dips JPfaz Trigo et aLll2009l L The presence of two relatively 
strong quasi-periodic signals (~ 112 days and ~ 190 days) 
in the light curve may still have certain relation with the or¬ 
bital period of NGC 5408 X-1. Such a feature can be caused 
by variations of the accretion rate due to a precession of the 
accre tion disc induced by the tidal force (iKotze fc Charlea 
l2012ll. or instabilities of disc warping driven by irradiation 
( Charles et al.l I2OO8I : lOeilvie fc Dubua [20o3) . Last but not 
least, the qua si-periodic signals may be associated with a 
precessing jet dFoster et al. l2of^ f. We therefore suggest the 
detected periods to be superorbital periods since these sig¬ 
nals are only prominent in a specific time interval of the 
long-term light curve although it is not entirely clear which 
mechanism induces these signals. According to previous in¬ 
vestigations between superorbital periodicities and orbital 


periods of Galactic X-ray binaries ( Kotze_^Charleg 2012l 'l . 


we support the conclusion given by 


Grise et al.l (l2013l L who 


report an orbital period of < 40 days. 


5.5 Concluding Remarks 

We have used the Swift archive to investigate the long¬ 
term X-ray variability of four ULXs. With timing analyses 
through the LSP, DPS and WWZ, we detect the following 
intriguing periodic signals embedded in their light curves: 
~ 380 days for ESO 243-49 HLX-1, and ~ 370 days and 
~ 110 days for M81 X-6. A brief summary of all potential 
periods is shown in Tabled 

Except for the main signal detected from ESO 243-49 
HLX-1, all other period candidates are only prominent in 
a specific time interval. The ~ 380-day periodic signal of 
ESO 243-49 HLX-1 corresponds to the separation between 
outbur sts, and it was clai med to represent an orbital modu¬ 
lation jLasota et ahllioill L Among the four ULXs discussed 
in this paper, only ESO 243-49 HLX-1 has a data set char¬ 
acterized by regular monitoring, which is adequate to be ex¬ 
amined with the HHT. This is also the first time that HHT 
has been applied to study long-term X-ray variability of a 
ULX. As compared with the DPS, the WWZ and the HHT 
are independent of the window time, and the Hilbert spec¬ 
trum provides a better resolution to discriminate the evolu¬ 
tion of each signal. We find that the separation between each 
outburst varies over time, and suggests this main signal to 
be a superorbital period. This may also suppor t the idea of 
a precessing X-ray beam (iKing fc Lasotall2014l) although a 
chan ging orbital period could also be possible JGodet et al.l 
I2OI4II . At the time of writing, the new outburst has started 
in 201 5 early-January, mo re than 460 days since the last out¬ 
burst (iKong et al.ll2013 L More monitoring observations in 


© 0000 RAS, MNRAS OOO.fTHTSl 






























































14 L. C.-C. Lin, C.-P. Hu, A. K. H. Kong, D. C. C. Yen, J. Takata and Y. Chou 


Table 2. Summary of our detected signals significant at > 99% 
against white and red noise 


Source 

Period/MJD range 

Possible Origin 

ESO 243-49 

HLX-1 

~ 380 days 
(--MJD 55050-56561) 

orbital period? 
superorbital period? 

Holmberg IX 
X-1 

625 days* ? 

(-^MJD 54803.7-57439.4) 

superorbital period ? 

M 8 I X-6 

~ 370 days 

(~MJD 54924.6-56554.7) 

~ 110 days 

(~MJD 55889.6-56554.7) 

superorbital period ? 

superorbital period ? 

NGC 5408 X-1 

None 

— 


* Too long to be definitely confirmed with current investigations 


the future are required to distinguish among different mod¬ 
els. 

The period candidates of other three ULXs may range 
from ~ 100 days to ~ 600 days. Apart from noise and 
artifacts, all the candidate periods are only significant in 
a specific epoch. This suggests that they are not associ¬ 
ated with any stable mechanism such as orbital motion. 
Instead, such long-term (> 100 days) X-ray quasi-periodic 
variations are likely related to superorbital periods that 
are thought t o be due to radiation- driven warping of ac¬ 
cretion discs llOgilvie fc Dublin 200lh or tidal in teraction- 
induced disc precession ilWhitehurst fc Kinelfl991^ . Alterna¬ 
tively, mass transfer rate-related events such as X-ray state 
changes and disc instability can also cause long-term modu¬ 
lations llKotze fc Charle^l2012l l. In particular, there are two 
intermittent quasi-periodicities for both NGC5408 X-1 and 
M81 X-6, suggesting that the quasi-periods are changing 
or evolving. They resemble some Galactic X-ray binaries 
th at show similar behavi our (e.g., Gyg X-2 and SMG X- 
1; iKotze fc Charle3 bOld l and it has been suggested that 
a warped disc could lead to an unstable steadily precess- 
ing d isc, causing quasi-periodic behaviour llOgilvie fc Dubua 
I2OOII I. We note that there are many uncertainties on the 
physical parameters of ULXs. To determine the origin of su¬ 
perorbital periods of ULXs, one has to know at least the 
mass ratio between the companion and the compact star 
{q = Mc/Mx) and the binary separation. Unfortunately, 
it is very difficult to get these parameters for ULXs. For 
the three ULXs discussed here (i.e., excluding ESO 243- 
49 HLX-1), only M81 X-6 has better constraints on the 
black hole mass and the nature of the companion. The 
masses of the black hole and companion star are estimated 
{Mx = 18 Mq, Me = 23 Mq) such that q can be de¬ 
rived. In this case, we can rule out a tidal interaction- 
induced disc precession sce nario that requires q < 0.25 — 0.33 
l| Whitehurst fc Kind Il99ll l. For a warped disc, the binary 
separation and the mass ratio suggest that M81 X-6 lies in 
the intermediate instability zone fo r radiation-driven warp - 
ing in X-ray binaries (see Figure 1 of iKotze fc Gharl^l2012ll . 
The quasi-periodic variability may represent the switching 
timescale between a warped disc and a flat disc. 

If all of these intriguing signals do represent superor¬ 
bital periods, the real orbital perio dicities of these system s 
are expected to be much shorter dKotze &: Gharle^ I 2 OI 2 II . 


This will give us constraints on the binary separation as well 
as the black hole mass. In order to understand the physical 
nature of the quasi-periodic modulations of the four ULXs 
discussed here, a regular long-term X-ray monitoring is re¬ 
quired in the future. More importantly, we demonstrate that 
dynamic timing analysis is important to investigate the long¬ 
term X-ray behaviour of ULXs especially when the periodic 
behaviours are intermittent or varying. 
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